7-Mn at the border between weak and strong correlations 
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We investigate the role of magnetic fluctuations in the spectral properties of paramagnetic 7- 
Mn. Two methods are employed. The Local Density Approximation plus Dynamical Mean-Field 
Theory together with the numerically exact quantum Monte-Carlo solver is used as a reference for 
the spectral properties. Then the same scheme is used with the computationally less demanding 
perturbative spin-polarized fluctuation-exchange solver in combination with the Disordered Local 
Moment approach, and photoemission spectra are calculated within the one-step model. It is shown 
that the formation of local magnetic moments in 7-Mn is very sensitive to the value of Hund's 
exchange parameter. Comparison with the experimental photoemission spectra demonstrates that 
7-Mn is a strongly correlated system, with the Hubbard band formation, which cannot be described 
by the perturbative approach. However, minor change of parameters would transform it into a 
weakly correlated system. 



PACS numbers: 71.20.Be, 75.10.Lp, 79.60.-i 

I. INTRODUCTION 

Itinerant electron magnetism has been a long-standing 
problem in solid state physics, due to the coexistence of 
itinerant and localized features in the behavior of mag- 
netic d-electronsii^i^. The standard ab initio approach to 
many-electron systems, i.e. the Density-Functional The- 
ory (DFT) in Local Density Approximation (LDA), or 
Generalized Gradient Approximation (GGA)^, describes 
the ground state properties of magnetic 3d metals rea- 
sonably well, but has serious problems with the descrip- 
tion of excitation spectra, e.g., photoemission spectra of 
nickel, and also with the electronic structure and mag- 
netism of iron, cobalt, and nickel at finite temperatures^. 
Correlation effects are essential also in other 3d itinerant- 
electron magnets, such as half- metallic ferromagnets^. 

A qualitative theory of itinerant-electron magnets tak- 
ing into account the duality of localized and itinerant 
behavior was developed in a framework of the Hubbard 
modelA. While this theory gave us a deep insight, it is 
insufficient to describe qualitatively properties of specific 
magnetic materials taking into account peculiarities of 
their real electronic structure. 

In the last decade much progress has been made with 
the introduction of the Dynamical-Mean Field Theory 
(DMFT)ii^, which allows for a quantitative solution of a 
lattice many-body problem with the inclusion of all lo- 
cal quantum dynamical fluctuations. The LDA+DMFT 
scheme, arising from the combination of the DMFT and 
standard DFT-LDA, has been successfully applied to de- 
scribe the electronic structure and magnetic properties 
of Fe and Ni^iSiiiS. In particular an excellent quantitative 
description of Ni was obtained (see also Refs. [Tllfl3 ). 



while for Fe the situation is slightly worse. This is not 
surprising, since in Fe the role of non-local effects (beyond 
DMFT) is more important than in Ni, as follows from di- 
rect photoemission data on the effective mass renormal- 
ization in different points of the Brillouin zone^^. 

Besides Fe and Ni, the electronic structure and mag- 
netic properties of 7-Mn are especially interesting. This 
phase appears between T — 1368 K and T — 1406 K, but 
can be stabilized down to room temperature through the 
addition of a small amount of impurities'^'^ or as layer- 
by-layer deposition on Cu3Au(100) ^^i^^ . Below the Neel 
temperature, about 540 K, 7-Mn becomes antiferromag- 
netic, which is accompanied by a tetragonal distortion 
into the fct structur o^^d"^ . Among all transition met- 
als, 7-Mn shows the largest discrepancy between the 
computational DFT-LDA (or GGA) predictions for the 
lattice constant and bulk modulus, and the experimen- 
tal dat a^^d^i^o . Tetragonal distortion and anomalies in 
the bulk modulus are good indications of strong elec- 
tron correlations, and Zein^^ has pointed out that the 
correlation effects in the 3d-shell are stronger at half- 
filling. In fact we have recently proved that inclusion of 
proper many-body effects through LDA-I-DMFT leads to 
correct equilibrium volume and bulk modulus^^. More- 
over, Biermann et alji^ have successfully applied the 
LDA+DMFT scheme to the angle-resolved photoemis- 
sion spectra (ARPES) of 7-Mn, emphasizing the forma- 
tion of a three-peak structure characteristic for strongly 
correlated systems. These calculations were carried 
out with the numerically exact quantum Monte-Carlo 
(QMC) solver, but, unfortunately, this technique is com- 
putationally quite expensive, especially, if one is inter- 
ested in the regime of relatively low temperatures, such as 
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room temperature. In addition, it uses a truncated multi- 
band Hubbard Hamiltonian, and involves many technical 
problems, e.g., with analytical continuation (for a more 
detailed discussion, see Ref. ;2i). 

In the present paper we address the question to which 
extent such a regime of strong correlations can be stud- 
ied by means of a simpler implementation. In partic- 
ular, the spin-polarized T-matrix fluctuation-exchange 
solver (SPTF)^i^ was implemented to treat the prob- 
lem of magnetic fluctuations in transition metals, and 
has been successfully applied to the ferromagnetic phases 
of Fe, Co, NiiiiSii^ and the anti-ferromagnetic phase of 
7-Mi*i^, as well as to half-metallic ferromagnet&§. It is 
quite stable, computationally cheap and deals with the 
complete four-indices interaction matrix. On the other 
hand, its perturbative character restricts its use to rela- 
tively weakly, or moderately, correlated systems. 

Not surprisingly, the SPTF performs well when start- 
ing from a spin-polarized solution, since the spin-splitting 
contains already the main part of the exchange and corre- 
lation effects. Conversely, the direct application of SPTF 
to a non-magnetic phase can create stability problems 
since we are trying, in this case, to attribute strong, and 
essentially mean- field, effect of formation of local mag- 
netic moment to dynamical fluctuations around non-spin- 
polarized state. This is appropriate when one uses the 
QMC method, which has no formal restrictions on the 
amplitude of fluctuations, but seems problematic for the 
perturbative approaches. 

As a way to weaken such a limitation we propose a 
combination of SPTF with the disordered local moment 
(DLM) approach^^'^^. As already shown for the case of 
actinides^^ the inclusion of the fluctuations of randomly 
oriented local moments can improve drastically the de- 
scription of energetics in the paramagnetic phase. One 
can hope therefore that it allows us to extend the range 
of apphcability of SPTF. 

As we show here, the results appear to be strongly 
dependent on the formation of local moments, which is 
driven by the Hund's exchange. For realistic values of 
the Coulomb interaction 7-Mn is shown to be a material 
right at the border between weak and strong correlations. 

II. METHOD 

In the LDA+DMFT scheme we correct the one-particle 
LDA Hamiltonian with an additional two-particle term, 
which explicitly describes the local Coulomb repulsions 
between a given set of correlated orbitals. The latter 
are usually associated to electrons in narrow bands, not 
properly described in the LDA. Naming R the lattice 
sites, the Hamiltonian becomes 

H = ^LDA + ^ X! X! ^1:2,3,4Cr,iC^,2CR,4CR,3, (1) 
R 1,2,3,4 

where the second sum runs over the quantum numbers 
of the chosen orbitals. The system in eq. (H]) repre- 



sents an "effective" Hubbard model, whose solution is 
a complicated many-body problem. The main observ- 
able is the local Green's function G(iw„), which can be 
obtained from the self-energy S(ia;„) by means of the 
Dyson equation. Here and below iujn are the Matsub- 
ara frequencies for finite temperature many-body formal- 
ism. In the DMFT the self-energy is obtained by map- 
ping the original lattice problem onto a problem of an 
atomic site coupled to an electronic bath, i.e. the effec- 
tive Anderson impurity model. The equality of the local 
Green's function in the original and the fictitious system 
assures the self-consistent condition that determines the 
bath. In the DMFT the self-energy S(ia;„) is purely lo- 
cal in real space, and consequently is k-independent in 
the momentum space, which becomes exact in the limit 
of infinite space dimensionality. Finally, notice that the 
Hamiltonian in eq. H]) is constructed on top of the LDA 
Hamiltonian, which should already contain all the effects 
of the Coulomb repulsions in the form of a mean-field. 
Thus, we remove from the self-energy those contribu- 
tions already included to avoid counting them twice. For 
treating metals the most common choice for the so-called 
" double counting" correction is the static part of the self- 
energy^^'^i, averaged over the orbital indices for each spin 
channel. 

In the present paper we have used two different 
LDA+DMFT implementations. The first one is the 
LDA+DMFT code of Refs. [Ulli, based on the full- 
potential linear muffin-tin orbital (FP-LMTO) method^-. 
We have used it for the QMC simulations, after having 
implemented the Hirsch-Fye Quantum Monte-Carlo algo- 
rithm of Refs. The second code is based of the full- 
potential Korringa-Kohn-Rostoker (FP-KKR) multiple- 
scattering theory, and has been used for the DLM+SPTF 
simulations^. The LDA+DMFT simulations with the 
SPTF solver were made using both these codes, and the 
results are practically identical, consistently with our re- 
cent findings for the ground-state properties of Ni and 
M»i^. 

The reason behind this twofold approach lays in the 
simplicity of the corresponding implementations. In 
QMC the interval [0, (3] on the axis of the imaginary 
time T is discretized in L time-slices. The fully inter- 
acting Green's function is evaluated only for the cho- 
sen T, and then transformed to Matsubara frequencies. 
Due to the statistical error, typical of Monte-Carlo meth- 
ods, and to the Trotter error, associated to the imag- 
inary time-discretization, the resulting G{iLUn) has bad 
analytical properties. While for FP-LMTO this is a mi- 
nor issue, in FP-KKR it is a serious problem, since the 
LDA+DMFT cycle involves analytical continuation from 
Matsubara axis to a semi-circular contour-Si. 

For DLM+SPTF the situation is inverted. In DLM 
the itinerant electrons form self-maintaining "local mo- 
ments" which are analogous - but physically different - 
to the localized spins of the Heisenberg model. Being 
the orientations of the moments at the sites i, we 
can describe the system through the generalized grand- 
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canonical potential rj({ei}). Then a mean-field approxi- 
mation of the true potential is constructed as expansion 
around a single-site spin Hamiltonian^^: 

f^o({e^}) = -^h,e,. (2) 

i 

The self-consistent parameters hi define a set of proba- 
bilities Pi{ei) of finding the moments oriented towards 
ei. Explicit calculations can now be made through stan- 
dard methods used for substitutionally disordered alloys, 
as the Coherent Potential Approximation (CPA), which 
has straightforward implementation in the multiple- 
scattering theory of KKR^^. On the other hand im- 
plementations of DLM into LMTO are more cumber- 
some, and involve the passage to a Green's function 
formalism^'^. 

Finally, for a direct comparison between theory and 
experiment, photoemission spectra are produced within 
the so-called one-step model2ii2£, which has been recently 
implemented in combination with the LDA+DMFT us- 
ing the KKR methodii. The main idea is to describe the 
excitation process, the movement of the electron towards 
the surface, and the final escape from the surface, as a 
single quantum-mechanically coherent process, which is 
comprehensive of all the multiple scattering events. 

To allow comparison with previous LDA-I-DMFT re- 
sults and experimental data, we assumed a face centered 
tetragonal structure with lattice constant a = 7.143 a. u. 
and a tetragonal distortion corresponding to c/a = 0.93. 
We have included 4s, Ap and 3d in the valence electrons; 
all the other states were considered as core states. In the 
FP-LMTO basis functions with three tail energies k were 
used, while in FP-KKR scattering matrix elements up to 
Imax — 3 were included. Only negligible differences could 
be found in comparing the LDA results. 

In the LDA4-DMFT simulation, the local Hubbard in- 
teraction was applied to the 3d electrons. In this case 
the ?7-matrix of eq. ([T]) can be parametrized by means of 
only two parameters U and J, which were varied among 
a wide range of reasonable values. In fact the determi- 
nation of J7 is a quite serious problem, and usually semi- 
empirical values are assumed^. We have studied config- 
urations corresponding to different parameters through 
QMC at temperature ranging from T ~ 2000 K down to 
T ~ 500 K. Convergence in the number of time-slices has 
been checked, and for high U values or low temperature 
a number of L = 128 time slices has been used. The 
number of Monte-Carlo sweeps has been set to 500000, 
while its convergency has been checked up to 10*". 

III. RESULTS 

In fig. [T] we show the density of states of the 3d elec- 
trons obtained through maximum-entropy method. In 
the bottom of the figure we have used an exchange pa- 
rameter J = 0.9 eV, and the results reproduce quite rea- 
sonably the three-peak structure of Ref. TgI. By turning 
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FIG. 1: (colour on-line) Density of states of the 3d electrons of 
7-Mn from QMC in the LDA-I-DMFT scheme for T = 2000 K. 
In the top plot results for J = 0.75 eV and U changing from 
3 eV to 6 eV are shown. In the bottom plot results for J = 
0.9 eV and same U as above are shown. Notice the big impact 
of small variations of J on the final spectrum. 



on the Coulomb repulsion [/, two split Hubbard bands 
are formed and a quasiparticle resonance appears close 
to the Fermi level. In comparison to the cited results, 
our central peak is almost completely smeared out be- 
cause of the high temperature, but still its foot-print 
can be seen in the spectrum. In the top of fig. [1] we 
have presented the density of states for the same simula- 
tions, but with a slightly different exchange J = 0.75 eV. 
Strong differences can be observed: the system can form 
Hubbard bands only for a very strong U and the en- 
ergy separation between them is smaller than the bare 
U. While for high J the electrons are observed to redis- 
tribute equally among all the orbitals in a way to mini- 
mize the QMC double occupation, for small J this does 
not happen and the occupations of the 3d orbitals stay 
closer to the original LDA picture. This behavior is re- 
lated with the formation of local magnetic moments, as 
it will be discussed below. Note that the transition from 
J = 0.9 eV to J = 0.75 eV looks small but, actually, it is 
very essential since it increases the difference C/ — 3 J by 
a factor of 2.5 which is the energy of the Coulomb inter- 
action between electrons with parallel spins in different 
orbital states^-. 

The problem is that the transition between these 
two regimes happens close to the real physical value of 
J. For all the 3d transition metals the semi-empirical 
atomic value J = 0.9 eV is commonly used^, but some 
calculations suggest that a lower value can be more 
adequate^ii^^. 

More information about the different physical situa- 
tions can be obtained by looking at the square of the 
local moment (M^), which in QMC is associated to the 
local spin correlation function. This moment depends 
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FIG. 2; (colour on-line) (a) Variation of the squared local 
moment (Mf ) as a function of T for chosen values of U and J. 
(b) Variation of the squared local moment (Af^) as a function 
of [/ at T = 2000 K. 



strongly on the temperature, as we can see m fig. [2l^a) 
for three pairs of parameters U and J. For very high 
temperatures T, i.e. for T bigger than J7, we expect the 
system to show a contribution to the effective Curie con- 
stant of about 0.5 for each orbital, which corresponds to 
local spins decoupled from each other—. In fact this is 
the tendency we observe in our calculations, not shown 
in the figure. For low temperature, i.e. in the interval 
between 2000 K and 500 K, the physics of manganese is 
strongly dependent on J. For low J the spin fluctuations 
are decreased because of the presence of strong orbital 
fluctuations^, which become important when the tem- 
perature becomes lower. The electrons tend to behave as 
a Pauli paramagnet, and the local moment decreases to- 
gether with the temperature. Conversely for high J the 
electrons tend to localize in each band independently, re- 
sulting in the suppression of the orbital fluctuations. As a 
result a strong fluctuating local moment can form. These 
different tendencies have been well studied in the two- 
orbital Hubbard model^*^ , and it is interesting to see how 
they can be found in a real material with flve (almost) 
degenerate orbitals and realistic hybridization. Notice 
that the experimental local magnetic moment is equal to 
2.3 Bohr magnetons. In fig. [IJb) we can see the trend 
of the local moments for different values of U in the two 
different regimes. 

In fig. [3] a comparison of all the significant density 
of states of our work is presented. The upper panel (a) 
and the middle panel (b) display respectively the bare 
LDA results and the LDA-hDMFT results with the QMC 
solver for J7 = 3eV and J = 0.75 eV or J = 0.9 eV. 
In the lower panel (c) the density of states obtained 
within the DLM approach is reported, both for LDA and 
LDA+DMFT with the SPTF solver ("DLM-I-SPTF"). 
While bare DLM describes the fluctuations in a very 
simple way, still a sort of three peak structure is ob- 
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FIG. 3: (colour on-line) Comparison of the significant density 
of states of the 3d electrons presented in the paper: (a) DFT- 
LDA (b) QMC with U = 3.0 eV and different values of J. (c) 
DLM and DLM-^SPTF for U = 3.0 eV and J = 0.8 eV. 



served. However the width of the Zd band is too big, 
since it descends from the single particle LDA density 
of states. In DLM+SPTF we can properly describe the 
shrinking of the Sd band, as already observed for Fe, Co 
and Ni2^. In comparison with the QMC calculations, 
the peak at around — 2eV is less pronounced. This de- 
pends on the perturbative nature of SPTF, which tends 
to shift the correlation effects related to the formation of 
non-coherent satellites. For example, the famous — 6eV 
satellite of Ni is positioned at about —8 eV. For Mn part 
of the spectral weight is transferred to the region be- 
tween —4 and — 6eV. While DLM+SPTF can reproduce 
very well the density of states for moderately correlated 
regimes, the appearance of the Hubbard bands observed 
in the QMC simulation for high J cannot be reproduced 
in this approach. 

On the basis of these considerations, we can look at the 
photoemission spectra reported in fig. [4j together with 
the experimental data of Ref. [l^. The spectra are con- 
sidered at normal emission along the (100) direction, and 
for photon energies ranging from 14 eV and 70 eV, which 
corresponds to the V — X path in the Brillouin zone. In 
the experimental data we can observe a pronounced peak 
close to the Fermi energy and a non-dispersive "Hub- 
bard band" at around — 2.8eV, whose intensity increases 
for high photon energies. The DLM simulation can de- 
scribe these features only partially. First the "quasipar- 
ticle peak" close to the Fermi energy seems to be much 
broader than its experimental counterpart. In second 
place a weak Hubbard band appears at around — 4.5eV, 
which is an energy significantly higher than — 2.8eV. In 
addition at high photon energies the signal is mixed with 
a stronger contribution coming from the sp band, which 
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FIG. 4: (colour on-line) Photoemission spectrum along the (100) direction in normal emission for photon energies between 14 eV 
and 70 eV, which corresponds to the F — X path in the Brillouin zone. The experimental data of Ref. [l^ are compared with 
calculations within the one-step model based on DLM and DLM+SPTF. In the spectra of the column labelled DLM+SPTF II 
the imaginary part of the self-energy coming from DMFT has been removed. 



makes it difRcult to identify this feature with an effective 
Hubbard band. The use of the DLM+SPTF approach 
improves essentially the description of the states near the 
Fermi energy but smears completely the Hubbard band. 
Probably, this is due to an overestimation of the imag- 
inary part of self-energy for larger excitation energies. 
Indeed, if one takes into account only the real part of the 
self energy (DLM+SPTF II) the description of the pho- 
toemission spectra can be drastically improved. This re- 
minds of a situation with the GW calculations where ne- 
glecting the electron state damping sometimes essentially 
improves the results^. Now we clearly observe a non- 
dispersive feature at —3.5 eV, which is separated from sp 
contributions, and can be surely identified as an effect 
of strong correlations. Nevertheless the peak appears at 
higher binding energy than in the experiment, due to 
the perturbative nature of SPTF. Finally notice that our 
simulations emphasize that the peak at the Fermi level 
possesses a strong surface contribution, visible in all the 
spectra as a small shoulder, and particularly sensitive to 
the modeling of the surface barrier of the photoemission 
process. 

In conclusion real manganese is a strongly correlated 
metal, with Hubbard bands^. In such cases, instead 
of the cumbersome and computationally expensive QMC 
calculations, the DLM+SPTF approach can be used to 
describe adequately the energy spectrum close enough to 



the Fermi energy, but not in the whole energy range. For 
moderately correlated systems, which were modelled here 
by variation of J, the DLM+SPTF approach turns out 
to be adequate for the whole spectrum. 
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